Stochastic resonance in a single electron turnstile 
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In this paper, we report stochastic resonance (SR) in a single electron turnstile. It has been 
known that SR emerges by the cooperation of a weak periodic signal and noise in a bistable system. 
A periodic signal produces switching between two stable states by lowering the potential barrier. 
Even if the amplitude of the signal is not large enough to make the potential barrier disappear, 
the switching can occur with the help of noise. As a result, an output signal-to-noise ratio (SNR), 
which shows how much the switching is synchronized with the signal, can be enhanced by increasing 
the noise. We have formulated the SR in a single electron turnstile into the adiabatic two-state 
theory of McNamara-Wiesenfeld and confirmed its manifestation by numerical simulation. These 
results could be applied to detection of the periodic motion of very small electric charge under noisy 
environment. 



I. INTRODUCTION 

Stochastic resonance (SR) is a nonlinear phenomenon 
whereby the addition of noise can enhance the detection 
of weak stimuli. An optimal amount of added noise re- 
sults in the maximum enhancement, whereas further in- 
creases in the noise intensity only degrade detectability. 
Although a large number of phenomena ranging from 
physics and engineering to biology and medicine have 
been studied (see the good review [T]), the essential ingre- 
dients for SR consist of a threshold, subthreshold stim- 
ulus and noise in nonlinear systems. In this study, we 
especially consider a bistable system with a weak peri- 
odic signal and random noise. Experimental observation 
of SR in such bistable systems involved a Schmitt trigger 
circuit and a radio frequency superconducting quan- 
tum interference device (rf-SQUID) li3^ [4 . An rf-SQUID 
utilizing SR could improve its performance by reducing 
its environmental noise sensitivity. It may be naturally 
expected that SR in a single electron circuit can be real- 
ized since the duality between an rf-SQUID and a single 
electron box (SEB) has been confirmed j5j. The SR in 
an SEB must be useful for detecting small charges un- 
der noisy environment. An rf-SQUID is composed of a 
Josephson junction, a superconducting loop inductance, 
and a bias current source in parallel. The number of flux 
quantums nf trapped in the loop is controlled by chang- 
ing the bias current. On the other hand an SEB is com- 
posed of a tunnel junction, a gate capacitance and a bias 
voltage source in series. The number of excess electrons 
He trapped in an island, which is the electrically isolated 
region between the tunnel junction and the gate capaci- 
tance, is controlled by changing the bias voltage. A per- 
fect duality seems to exist, however, it is not true about 
bistability. While an rf-SQUID has bistability related to 
Uf due to hysteresis property of a Josephson junction, 
an SEB does not have bistability related to ng because 
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one tunnel junction does not show hysteresis property. 
Then let us consider a single electron turnstile j6j which 
is given bistability with four tunnel junctions in series. 
In the following sections, we formulate the SR in a single 
electron turnstile into the adiabatic two-state theory of 
McNamara-Wiesenfeld [2| and confirm its manifestation 
by numerical simulation. 



II. THEORY 

A. Analysis of the dynamics 

Figure [T] shows the equivalent circuit of a single elec- 
tron turnstile and its stability diagram. The stability di- 
agram is a two-dimensional map of the stable states of a 
circuit jSj. The x- and y-axes represent two bias voltages, 
namely the normalized gate voltage and the normalized 
bias voltage, respectively. At the operating points sur- 
rounded by each diamond in Figjljb), the change of n, 
which signifies the number of excess electrons in the cen- 
tral island, leads to increase of free energy. Then elec- 
tron tunneling is inhibited (Coulomb blockade [5]) and 
the circuit remains stable states. Thus the gray regions 
in which two neighboring diamonds are overlapped are 
bistable regions: two different stable states with differ- 
ent n have been realized exclusively. If the gate voltage 
Vg{t) = Ay sm{2n fst) + Vgo takes the operating point cut 
across the bistable region periodically like A-^ B^ A 
■ ■ ■ in Fig TJc) , n varies ^ 1 — > — > • • • . Con- 
sequently such single electron transportation synchro- 
nized with the frequency fs can be realized and then 
current Is = efs flows. This property can be utilized 
for a current standard. Notice that the periodic mo- 
tion of an electric charge Qs{t) = AgSm{27T fst) polarized 
in the central island corresponds to the ac-component 
Vs{t) — Ay sm{2TT fst), where A^ ~ Aq/Cg and Cg is a 
gate capacitance. While the periodic motion of an out- 
side electric charge near the central island produces the 
synchronized current flow in the same manner even with- 
out the ac-component of Vg. Therefore a single electron 
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FIG. 1: A single electron turnstile (a) and its stability dia- 
gram at zero temperature (b). The stability diagram enlarged 
near the bias point (c). An alternating gate voltage Vg{t) 
swings the operating point. 



turnstile can be used as a detector of the periodic motion 
of electric charge. Let us consider Qs{t) and n{t) as in- 
put and output signals, respectively. Bistability related 
to n{t) leads to SR with the help of some external noises 
such as fluctuation of voltages, fluctuation of background 
charge motion. 

Let us suppose the circuit is biased at the point 
(VgO, VbO - e) as shown in Fig|TJc), where Vgo = e/{Cg + 
Cext), VbO = e/2(C + Cext), e is elementary charge, C is 
a tunnel capacitance, Cext = C{C/2 + Cg)/(3C/2 -I- Cg) 
and e is a small positive constant. The input charge 
signal Qs(t) forces the operating point of the circuit to 
fluctuate with the amplitude Aq/Cg and the frequency 
fs around the bias point. We assume that Gaussian 
charge noise QN{t) with cutoff frequency /c is added 
to Qs{t)- Then the input gate voltage is given as 

Vg{t) = VgO+{Qsit)+QN{t)}/Cg = VgO + Vs{t) + Vn (t) . 

Switching of n{t) occurs a,t Vg — Vto{n, = 0^1) and 
Vg = Vti{n=1^0), where 
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2Vb{C+Cext)/iCg + Cext), 



2Vb{C+Cext)/{Cg + Cext), 



(1) 



and Vb is a bias voltage. Note that the signal frequency 
fs has to meet the "adiabatic limit" l/fs ^ Tt, where 
Tt is the tunneling time of electrons, in order to get the 
circuit settle in equilibrium states n =0 or 1. If tunneling 
time is negligibly small, the change of excess electrons 
accompany with tunneling at zero temperature is given 
as follows, 

n{t+) ^ ^|s.gn(z) -I- l|, 

z = Vg{t+) - [n{t^)Vn + {I - n{t^)}Vto] 
= {Vto-Vti)n{t^)~Vto + Vg{t+), (2) 



where sgn(z) = 1 if z > and sgn(z) = — 1 if z < 0. 
Subscripts + and _ for t mean after and before tun- 
nehng, respectively. The behavior defined by Eq.Q is 
similar to that of a Schmitt trigger circuit. Hereafter 
we assume that operating temperature T is not zero and 
then switching occurs stochastically in the vicinity of Vfo 
and Vti- Hence n(t) should be treated as a continuous 
variable so that it indicates an expected value. Following 
the procedure derived for the Schmitt trigger circuit [7], 
the time evolution of n(t) decaying with a time constant 
Tt is given as follows. 
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n - -{tanh(/3z) -f 1} 



(3) 



where /3 is the gain parameter related to T . Integrating 
the right side of Eq.([3| produces the potential U{n,t) 
which governs the dynamics of n{t). 



U{n,t) 



Vti) 



X \ncos\y[P{{Vto-Vn)n-Vto + Vg{t)}] .(4) 



The signal Vs{t) = Qs{t)/Cg modulates U{n,t) period- 
ically and lowers the potential barrier between the two 
stable states n = and 1. Even if the amplitude of 
the signal is not large enough to make the potential bar- 
rier disappear (subthreshold signal) , switching can occur 
with the help of the noise Vat (t) ^ Qn (t) /Cg ; this is the 
source of SR. 



B. Stochastic resonance — adiabatic two-state 
theory 

Since the dynamics of a single electron turnstile is de- 
scribed with the double- well potential U{n,t), we can 
apply the adiabatic two-state theory 7] in which SR in 
bistable systems have been studied. An output signal-to- 
noise ratio (SNR) shows how much the switching between 
two stable states is synchronized with the input signal. 
The addition of noise degrades the SNR in linear systems, 
while added noise results in the enhancement of the SNR 
in nonlinear systems in which SR occurs. Therefore the 
SNR versus input noise profile is sometimes taken to be 
the hallmark of SR. The SNR is computed by integrating 
the power spectral density (PSD) over the peak centered 
at the signal frequency fs , and dividing by the mean noise 
power around fs and is found in dB [7], 



SNR= 10 log 



A^A 



G is the processing gain given as follows, 

1 [Ef 



G = 



FFT 



wf 



(5) 



(6) 



where Wi is the window coefficient multiplying the zth 
sampling in the time series and NppT is the number of 
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sampling points in fast Fourier transform (FFT) . A is the 
width of a frequency bin and is equal to fsampie/NppT, 
where f sample is the sampling frequency. The ratio S/N 
is written in the form [7] 



S_ 
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ao and ai are calculated as follows, 

ao = 2W^(?7 = 0), 
dW 

.-2-(,.0). 



(7) 



(8) 



W{t) = f{^i + Tj) is the transition rate out of the stable 
states, where /i is a parameter related to a potential bar- 
rier and r] = r]o cos(27r/si) is a stimulus modulating the 
potential periodically. 

According to the previous results of Tsironis 
et al. [To] [11] , the transition rate obtained for a double- 
well system in the presence of colored noise is written in 
the form. 
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where Dy is the variance of Gaussian noise VAr(i), r^r = 
l/(2/c) is the correlation time of VN{t), and tn has to 
meet the "strong color noise limit" rjv 3> Tf for the va- 
lidity of Eq.(|9]). V/vc(i) is the critical noise which makes 
the potential barrier disappear in the presence of the sig- 



nal Vs{t) as shown in Figl2l 

VN{t) = VNc{t) yields 
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Substituting Eq.([To]) for EqM yields W{n,Vs) (?7o cor- 
responds to Vg's amplitude and both and ai are 
calculated by Eq.([8| and written as follows. 
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Although we introduce U (n, t) using a continuous ap- 
proximation for n in average meaning, the possible value 
of n is discrete. Then the single electron turnstile system 
can not be mapped onto the double- well system exactly. 
Another different approach calculating W is derived from 
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FIG. 2: Bistable potential modulated by a subthreshold signal 
Vs and critical noise Vnc- Switching of n{t) can occur with 
the help of noise larger than critical noise. 



a Fokker-Planck equation [7] and ao and ai are calcu- 
lated as follows. 
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where Erf(u) is the error function, A is a parameter 
related to T and 7 is obtained in our case as 7 = 
A(Vto — Vti)/2. Note that ryo in Eq.(|7]) corresponds not 
to but to A^/{2DvY/'^- The integral in Eq.([l2]) can 



be calculated numerically. The strong color noise limit is 
also required for validity of Eq. ( 12 1 . 



III. NUMERICAL SIMULATION 

We have confirmed the SR effect obeying the theory 
mentioned above by numerical simulation based on a 
Monte Carlo method [12j. Parameters are chosen as fol- 
lows; C = 1.0 [aF], the tunnel resistances Rt — 100 [kfi], 
Cg = 0.5 [aF], VgO = 0.160218 [V], V& = 50 [mV], T = 
30 [mK], /, = 100 [MHz], = 125 [ps] (/, = 4 [GHz]). 
In our parameters, Tt is several tens pico seconds on aver- 
age. The PSD related to the time series of central island 
voltage is computed by using a 2048 point FFT with 
2 [GHz] sampling, applying a Hanning window to each 
segment (one segment is composed of 100-cycle waves). 
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by Eq.(ll) agree with the numerical results only in the 
regions wherein the SNRs increase. While the results 
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FIG. 3: SNR versus noise profile. Point, solid line and dashed 
line plots correspond to the simulation results, the theoretical 
results by Eq.(lll and Eq.(12l, respectively. /? = 4900 and 
A = 1.63 are determined by a method of least squares. 



and averaging 100 segments. For N in Eq.([5l, we use 
an interpolated average of the PSD from neighboring 10 
frequency bins. Averaging 100 ensembles for each noise 
variance Dy , we have obtained the SNR versus noise pro- 
file numerically. 

Figure |3] shows three results obtained by numerical 
simulation (point plots) , by Eq. (11) (solid line plots) and 
by Eq.([T2| (dashed line plots) for various input signal 
amplitude Aq = CgA^. (3 = 4900 and A = 1.63 are 
treated as free parameters and determined by a method 
of least squares because it is difficult to determine them 
analytically. SR effect, wherein the SNR passes through 
a maximum at certain noise strength, has been con- 
firmed from these results. Theoretical results obtained 



obtained by Eq.(12l can predict correctly the locations 



of the SNR maxima and the decreasing rate of the SNRs 
through the maxima. The conspicuous differences are the 
rise-up points and the magnitude of the SNRs especially 
with large Aq. This is because that the adiabatic two- 
state theory is valid on condition that signal amplitude 
is much smaller than a potential barrier. Another possi- 
ble reason for the differences is a strong color noise limit 
violation because can sometimes be comparable to r^v 
due to ffuctuation of the tunneling time of electrons. 



IV. CONCLUSION 

We have studied stochastic resonance in a single elec- 
tron turnstile both theoretically and numerically. These 
results could be applied to detection of the periodic mo- 
tion of very small electric charge under noisy environ- 
ment. One possible application is qubit detection in 
quantum computers fe.o.|13j. |14j ) instead of using an 
SET (Single Electron Transistor) [3] or an rf-SET (radio 
frequency SET) [T^. Even though it detects not quan- 
tity but motion of electric charge, the superior figure, i.e. 
noise tolerance, is highly attractive for detecting weak 
signals. 
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